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ABSTRACT 



We are conducting a planet search survey with HARPS since seven years. The volume-limited stellar sample includes all F2 to 
M0 main-sequence stars within 57.5 pc, where extrasolar planetary signatures are systematically searched for with the radial-velocity 
technics. In this paper, we report the discovery of new substellar companions of seven main-sequence stars and one giant star, detected 
through multiple Doppler measurements with the instrument HARPS installed on the ESO 3.6m telescope, La Silla, Chile. These 
extrasolar planets orbit the stars HD 1690, HD 25171, HD33473A, HD 89839, HD113538, HD 167677 , and HD217786. The 
already-published giant planet around HD 72659 is also analysed here, and its elements are better determined by the addition of 
HARPS and Keck data. The other discoveries are giant planets in distant orbits, ranging from 0.3 to 29 Mj up in mass and between 
0.7 and 10 years in orbital period. The low metallicity of most of these new planet-hosting stars reinforces the current trend for 
long-distance planets around metal-poor stars. 

Long-term radial-velocity surveys allow probing the outskirts of extrasolar planetary systems, although confidence in the solution 
may be low until more than one orbital period is fully covered by the observations. For many systems discussed in this paper, 
longer baselines are necessary to refine the radial-velocity fit and derive planetary parameters. The radial-velocity time series of stars 
BD -114672 and HIP 2 1934 are also analysed and their behaviour interpreted in terms of the activity cycle of the star, rather than 



> 

o 

m 
oo 
m 



£NJ ■ long-period planetary companions. 

Key words, stars: individual: HD 1690 , BD -1 14672 , HD 25171 , HD33473A, HD 72659 , HD 89839 , HD 113538 , HD 167677 , 
HD 217786 , HIP 21934 - stars: planetary systems - techniques: radial velocities - techniques: spectroscopic 



1 . Introduction self-luminosity than a planet of similar mass in an evolved sys- 

jj] ■ tem. The combination of both methods will permit calibrating 

C$ ■ The complementarity between planet search techniques will al- the atmosph eric models, i.e., the relationship between the mass 

low us in a mid-future to draw the general picture of extrasolar of the planet and its spec tral energy distribution (e.g., Baraffe 

systems. The transit method is strongly biased towards short- et a j (2003)' Fortney et al (2008)) 
period planets, while direct imaging will be limited to exploring 

the outer regions of close-by systems. Intermediately, the radial- In this context, the volume-limited sample conducted with 
velocity surveys, conducted for more than 15 years, are contin- HARPS since 2003 is one of the main sources of favourable tar- 
uously discovering planets in the wide period range between 1 ge ts for future direct imaging surveys. With a precision of about 
day and 15 years. In some more years, we will be able to study 2 m/s over 850 stars, it has allowed the discovery of 32 plan- 
the same extrasolar systems by means of the gravitational pertur- e ts so far (Pepe et al. 2004; Moutou et al. 2005; Lo Curto et al. 
bation inferred by the planets on their star (by the radial-velocity 2006; Naef et al. 2007; Moutou et al. 2009; Lo Curto et al. 2010; 
and/or the astrometric technics), and direct imaging with the aid Segransan et al. 2010; Naef et al. 2010). In this paper, we present 
of the next-generation high-contrast instruments. This will be the HARPS data collected for eleven stars in this stellar sample, 
possible only for a few nearby systems that contain outer gi- which gathers F2 to M0 main-sequence stars closer than 57.5 
ant planets, or for young systems when the planet has a higher pc from the Sun. Their radial-velocity variations is analysed and 
interpreted as the signature of planetary companions, except in 



Send offprint requests to: C. Moutou two cases where we show strong evidence for the detection of the 

* Based on observations made with the HARPS instrument on the stellar long-term magnetic cycle. In section 2, the parent stars are 

ESO 3.6 m telescope at La Silla Observatory under programme IDs described, while radial-velocity observations and their analysis 

072.C-0488(E) and 085.C-0019. are presented in section 3; a discussion follows in section 4. 



2 



C. Moutou et al.: The HARPS search for southern extra-solar planets 



2. Characteristics of the host stars 

The spectroscopic analysis of the stars has been conducted on 
the combined HARPS spectra, following the method described 
in Sousa et al. (2008). The method is based on the equivalent 
widths of Fe I and Fe II weak lines, by imposing excitation and 
ionization equilibrium assuming LTE. For this task the 2002 ver- 
sion of the code MOOG was used (Sneden 1973) and a grid 
of Kurucz Atlas 9 plane-parallel model atmospheres (Kurucz 
1993). Equivalent widths of spectral lines are estimated with the 
ARES code (Sousa et al. 2007), from which the values of T e ff, 
logg and [Fe/H] are derived. Actual errors on these parameters 
are obtained by quadratically adding 60 K, 0.1 and 0.04 dex re- 
spectively to the internal errors on T e fj, logg and [Fe/H]. These 
values were estimated by considering the typical dispersion in 
each parameter comparison plot with the infra-red flux method 
(Casagrande et al. 2006), presented in Figures 3 to 5 of Sousa 
et al. (2008). The stellar luminosity was calculated based on 
Hipparcos photometry (ESA 1997) and revised Hipparcos paral- 
laxes (van Leeuwen 2007), including the bolometric correction 
as calibrated in Flower (1996). Theoretical evolutionary tracks 
are compared with the measured values to estimate the mass, ra- 
dius and age of the stars with a Bayesian analysis 1 (da Silva et al. 
2006). The width and contrast of the cross-correlation function 
are measured to estimate the projected rotational velocity of the 
star. Finally, the activity indicator logR' HK is measured on indi- 
vidual spectra when the flux is sufficient at 390 nm (Santos et al. 
2002). 

Tables 1 and 2 summarize the measured and derived param- 
eters of the stars. For the reddest stars, the relationships between 
the cross-correlation function and the stellar parameters (partic- 
ularly vsim) are not valid (Santos et al. 2000, 2002), and those 
values are thus missing. 

3. Radial velocity data and proposed solutions 

Radial-velocity time series of our sample stars have been ob- 
tained with the spectrograph HARPS (Mayor et al. 2003), 
mounted at the 3.60m telescope in La Silla, ESO, Chile, mostly 
in the framework of the Guaranteed Time Observation pro- 
gramme (programme ID 072.C-0488), with additions from the 
follow-up volume-limited search programme 085.C-0019. 

We used the HAM mode with a spectral resolving power of 
1 15000, without simultaneous calibration of the velocity by the 
ThAr calibration lamp: the spectrograph is indeed stable enough 
for the precision of about 2 m/s targeted in this survey. Such a 
precision corresponds to a signal-to-noise ratio of 40 at 550 nm 
and to exposure times ranging from 2 to 15 min according to the 
stellar magnitude and external conditions, with a mean value of 
8 minutes. 

Each individual spectrum was then cross-correlated with a 
numerical mask, built from the stellar template better matching 
the stellar type, from M0 to F0 (Baranne et al. 1996). In addition 
to the position of the cross-correlation peak, the slope of the bi- 
sector line is measured on the stellar profile, indicating the level 
of photospheric activity that may affect the measurement of the 
peak position. The variation of the bisector span with the radial 
velocity is shown in Figure 13 for all stars of the presented sam- 
ple. Let us recall that an anti-correlation of the bisector span with 
respect to the peak position may depict a spot-related variation 
of the radial velocity, especially when the amplitudes of both 
variations are similar. Bisector variations may also be used to di- 
agnose long term radial-velocity signals induced by variations in 
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Fig. 1. Left: Radial- velocity measurements of HD 1690 obtained with 
HARPS against time. The orbital fit is overplotted as a black line. 
The residuals to the fit are given below. Right: the phase-folded radial- 
velocity curve is shown. The best fit to the data is shown as a plain line. 



the stellar magnetic cycle (Santos et al. 2010a), as well as varia- 
tions in the width of the cross-correlation function and temporal 
evolution of the activity tracer logR' HK (see section 3.3) . 

In the time series of stellar radial velocities presented below, 
a significant variation was detected over periods usually longer 
than several hundreds days. After checking for activity indica- 
tors, we present an analysis of the results in terms of the plane- 
tary orbital signature except for two cases where the long-term 
activity cycle seems the dominant factor. In a few cases, several 
solutions are equivalent and more data over monthes to years are 
needed to conclude. Radial-velocity time series and the proposed 
solution are presented in Figures 1 to 12. Individual measure- 
ments are available in electronic form at the CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u- 
strasbg.fr/cgi-bin/qcat?J/A+A/. Error bars for the derived pa- 
rameters were estimated from a Monte-Carlo bootstrap analysis 
of the data sets using the routines of yorbit (written by one of 
us, D.S.). 

3.1. Complete orbits 
HD1690 

Fourty-two HARPS measurements of the giant star 
HD 1690 were gathered from November 2003 to October 
2010, as shown in Figure 1. The presence of a few giant stars 
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HD25171 HARPS 
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Fig. 2. Radial- velocity measurements of HD 25 171 obtained with 
HARPS against time and the residuals to the best-fit model. 
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Fig. 3. Variations of x 1 as a function of various related parameters of the 
fit for HD 25171, using Markov-chain simulations. (Left) Correlation 
plot between the orbital period and radial-velocity semi-amplitude. 
(Right) Correlation plot between eccentricity and orbital period. The 
lines show 1,2 and 3-cr iso-contour plots. 



When the variation of the stellar radial velocity of 
HD 1690 is interpreted in terms of a Keplerian orbit of a planet, 
the best solution corresponds to a 532 day period. The resid- 
ual scatter of the fit is rms = 25.3 m/s. Since the parent star is 
a giant, one expects the residual scatter to originate from high- 
frequency stellar oscillations, with an amplitude that is standard 
for red giants (Sato et al. 2005). No significant period is found 
in the residual of the fit, which is also compatible with the com- 
mon behaviour of radial-velocity oscillation signatures in giants 
(Hatzes et al. 2007). In order to better reflect the actual velocity 
jitter and to get a reduced^ 2 closer to 1, we quadratically added 
25 m/s to the individual errors. This jitter amplitude is common 
for a K giant with B - V=1.35, a range where it may vary from 
20 to 100 m/s (Frink et al. 2001). The residuals show a slight 
negative slope, indicating a possible distant companion. The so- 
lution implying a linear drift reproduces the data better than the 
single planet with a probability of 75% (F-test of Pourbaix & 
Arenou (2001)). The orbit of the planet shows large eccentric- 
ity, e - 0.64, leading to a minimum mass of the companion of 
about 6 Mj up . Figure 1 shows the phase-folded data points and 
best-fit orbit. It must be noted that the solution shown in this 
analysis is not very robust: the removal of -12 data points sig- 
nificantly changes the period of the best-fit. The reason for this 
are several peaks of equivalent power in the periodogram. The 
presented solution is found by performing 10000 Monte Carlo 
simulations and including all data points. The complete solution 
is given in Table 3. It is possible that further observations modify 
the solution proposed here, or increase the number of planets in 
the system. 



HD25171 

The non-active star HD 25 171 was observed 24 times with 
HARPS between November 2003 and March 2010 with average 
errors of 3 m/s. A long-period variation is detected and attributed 
to the presence of a distant Jupiter-like planet. The standard devi- 
ation of the velocity drops by a factor 3 when a best-fit Keplerian 
orbit is removed from the data. A period of 5 years is found, cor- 
responding to a planet minimum mass of 0.95 Mj up (Figure 2). 
Since the semi-amplitude is low and the long orbit only partially 
covered by observations, we performed a Monte-Carlo Markov- 
chain analysis to support the detection. The most relevant out- 
put plots are shown on Figure 3: the x 1 distribution of the fit, 
against related parameters such as the semi-amplitude, the ec- 
centricity, and the orbital period. The correlation plots show that 
several equivalent solutions exist, but in a limited range of these 
parameters The complete solution is given in Table 3 from the 
bootstrap analysis that gives more conservative errors than the 
MCMC analysis. The scatter of the residuals to this fit is 3.4 
m/s, which is of the order of individual measurements. 



HD 72659 



in the HARPS volume-limited sample is due to a revision of 
the parallax from the original Hipparcos output catalogue (ESA 
1997) to the new estimates of parallaxes from van Leeuwen 
(2007). For star HD 1690, the parallax was revised from 43.42 
to 3.22 mas. 

The average formal accuracy of individual data points is 1.8 
m/s. Variations of the radial velocity have a standard deviation of 
1 10 m/s, while the bisector slope varies with a standard deviation 
of less than 20 m/s. 



A total of 36 spectra of HD 72659 have been gathered between 
February 2004 and June 2010, with an average error of 2.3 m/s 
(Figure 4, red points). The radial-velocity time series shows vari- 
ations with a maximum amplitude of 70 m/s, whereas the bi- 
sector slope is stable to 6 m/s. The full period is seemingly not 
covered by our observations. The low activity of the star and 
the absence of a correlation of the bisector span or activity trac- 
ers with the radial velocity exclude the possibility that activity 
alone explains the velocity variation. 
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A planetary companion has already been reported by Butler 
et al. (2003) in this system, with poor constraint on the orbit. A 
refined solution is published in Butler et al. (2006). The combi- 
nation of both Keck/HIRES and HARPS data reinforces the so- 
lution with more than 12 year span and in particular, the orbital 
period is now well determined. The offset between the Keck and 
HARPS data sets is left as a free parameter to the fit. Our so- 
lution is compatible with the one given in Butler et al. (2006), 
but with an error on the period seven times smaller. The best-fit 
solution shows the presence of a single orbiting 3.1-My„ p planet 
in a slightly eccentric orbit of period 3658 days (Figure 4). The 
residuals of the fit have a standard deviation of 3.7 m/s. No long- 
term trend or significant peak in the periodogram of residuals 
is found to explain this low-amplitude extra jitter. The complete 
solution is given in Table 3. 



HD72659 HARPS Keck 
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Fig. 4. Radial-velocity measurements of HD 72659 obtained with 
HARPS (red points) and Keck (blue points) against time, together with 
the best-fit solution of the combined data sets and the residuals. 



HD 113538 

We secured 29 HARPS measurements of star 
HD 11 3538 between February 2004 and June 2010; they 
have average errors of 1.8 m/s. The radial velocity variations 
have a peak-to-peak amplitude of 50 m/s over the 2268 day 
time span of the observations, and they do not correlate with the 
bisector span variations. The best-fit solution involves a single 
planet with a long period and very eccentric orbit, but residuals 
to the fit are large with regard to individual errors (standard 
deviation of 5.4 m/s), and with a poorly defined semi-amplitude. 
A better solution is depicted with a two-Keplerian orbit model, 
involving a Jupiter-like planet in a slightly eccentric orbit of 
1657 day period, plus a Saturn-like planet in a shorter eccentric 
orbit of 263 days. Then the residuals to the combined fit drop to 
1.8 m/s. The complete solution is given in Table 4 and Figure 
5. Simulations were performed to confirm that the two-planet 
solution is statistically more significant than the single planet 
+ long-term trend variation: the data are permuted around the 
one-planet solution, and the residuals are statistically compared 
with the two-planet model residuals. The probability that the 
two-planet model suits the data better is obtained by calculating 
the number of realisations where the x 1 obtained with the 



permuted data is smaller than with real data. The use of the 
genetic algorithm in this iterative process allows more reliability 
than false alarm probability tests on the periodogram. 

As an additionnal test and because the star is significantly 
active, we used the activity indicators to check whether the sig- 
nal is caused by activity. The FWHM of the cross-correlation 
function shows a moderate correlation to the Mount Wilson S 
activity index, which is based on calcium emission, with a cor- 
relation coefficient of 60%. Our previous analysis shows that the 
two-planet model is statistically robust, but we remain cautious 
because stellar activity probably affects the solution. 

HD 21 7786 

Fifteen spectra of HD 2 17786 were acquired from July 2006 to 
October 2009. They show very large velocity variations, with a 
peak-to-peak amplitude of 390 m/s when individual errors are 
2.1 m/s on average. The bisector span shows small variations 
that are not correlated with the radial velocity, the FWHM and 
the calcium activity index are not correlated, and the star is not 
active with logR' HK of -4.86. The velocity variation is thus not 
caused by activity. A Keplerian solution with a 1320 day period, 
0.4 eccentricity and 260 m/s semi-amplitude is the best-fit solu- 
tion to the data set (Figure 6). The residuals to the fit have an rms 
of 2.5 m/s. In order to further increase the significance of the fit, 
additional observations around periastron and apastron would be 
necessary. The minimum mass of the orbiting planet is 13 Mj up , 
within the overlapping region between giant planets and brown 
dwarfs. The complete solution is given in Table 3. 

3.2. Incomplete orbits 
HD 33473 A 

Twenty-six HARPS spectra of HD 33473 A were collected from 
October 2003 to March 2010 with an average error of 1.9 m/s. 
Large amplitude variations of the velocity have a standard de- 
viation of 200 m/s. A velocity variation caused by stellar ac- 
tivity is rejected by our data, because the bisector span is con- 
stant. The 2300-day time span covered by our observations is 
not long enough to distinguish between two solutions: an un- 
complete, single orbit with an ill-constrained period (16 to 38 
years), or a ~2800 day period orbit with a longer-term drift. Both 
solutions give a similar^ 2 . In another year, with continuing ob- 
servations of this system, we should be able to better constrain 
the system, because then the shorter period will be entirely ob- 
served. With such a difference in the derived period, the mass of 
the companion is evidently very uncertain, ranging from 6 to 29 
M j up . However, HD 33473 A has a known visual companion, HD 
33473B, with an angular separation of 10 arcsec and magnitude 
difference 3.4 in the optical. It is thus a K8-like star at 500 UA 
distance, at least, which corresponds to a period longer than 10 4 
years. The existence of this distant stellar companion strongly 
reinforces the justification of a planet model involving a long- 
term drift. This is the solution presented in this paper (Table 5 
and Figure 7). It gives an O - C residual scatter of about 2.4 m/s. 

HD 89839 

From February 2004 to June 2010, 39 measurements of the 
low-activity star HD 89839 were recorded. They show long-term 
radial-velocity variations with a peak-to-peak amplitude of 90 
m/s (Figure 8). No correlation is found between the bisector and 
the velocity, excluding the stellar activity as the origin of the 



C. Moutou et al.: The HARPS search for southern extra-solar planets 



5 



HD113538 



HARPS 



HD217786 



HARPS 



> 
CC 




-14 



3000 3200 3400 3600 3800 4000 4200 4400 4600 4800 5000 5200 5400 
JD - 2450000.0 [days] 



HARPS 



> 
cc 




- LlUlLlI - 

-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 I.I 




-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



Fig. 5. Radial-velocity measurements of HD 1 13538 obtained with 
HARPS against time, and for both planets, against their orbital phase. 



observed signal. Although the orbit is not regularly sampled by 
the observations, a period of 9.5 to 27 years days emerges, with 
some excentricity. As we are entering the descending section of 
the radial velocity curve, this solution will be confirmed and re- 
fined within the next four years. The current solution gives a 
semi-amplitude variation of 45 m/s, corresponding to a 3.9 M j up 
planet. The standard deviation of the residuals to the fit is 3.9 
m/s, similar to the average error bar of the measurements. The 
MCMC output distributions show the correlated behaviour of ec- 
centricity, semi-amplitude, and orbital period in Figure 9. 
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Fig. 6. Radial- velocity measurements of HD 217786 obtained with 
HARPS against time. 



HD 167677 

Twenty-six HARPS measurements of HD 167677 were collected 
between May 2005 and October 2010, with an average Icr error 
of 1.6 m/s. They show radial-velocity fluctuations at the level of 
15 m/s standard deviation, and with a periodicity of about 1800 
days. The activity indicator log R' HK is typical of a non-active star 
(-4.99). The bisector span and FWHM of the cross-correlation 
function do not vary in phase with the calcium index or the ve- 
locity. So the interpretation in terms of pure stellar activity is 
excluded. The best-fit Keplerian solution corresponds to a planet 
with minimum mass 1.36 Mj up in a slightly eccentric orbit, as 
shown in Figure 10. The residuals to the fit have a standard devi- 
ation of 2.6 m/s. The MCMC analysis allowed us to quantify the 
correlation between the orbital parameters, and Figure 1 1 shows 
the relation between eccentricity, period, and semi-amplitude. 
Additionnal data would be needed to get more precise param- 
eters of the orbit, and can be collected in the next years. The 
complete solution is given in Table 5. 

3.3. Long-period planets or magnetic cycles? 

3.4. BD -11 4672 

Nineteen HARPS measurements were obtained of star BD - 
1 14672 between July 2005 and June 2010, with an average error 
of 1.9 m/s (Figure 12). A peak-to-peak velocity variation of 28 
m/s is observed. The bisector slope shows a significant degree of 
anti-correlation (Fig. 13), and a visible reversal in the core of the 
calcium lines confirms that the star is active. However, the main 
radial-velocity variation has a significant peak of periodicity at 
1692 days, much too long for a rotation period, and no signifi- 
cant peak in the range 5-50 days. An attempt has been made to 
correct for the activity as in Melo et al. (2007) and Boisse et al. 
(2009) by linearly fitting the bisector-velocity correlation and 
subtracting this trend from the velocity measurements. The re- 
sulting solution is not significantly different. In addition, Figure 
14 shows a significant degree of correlation between the full- 
width-at-half maximum of the cross-correlation function, and 
the calcium activity index (Mount Wilson S value). The stel- 
lar lines therefore appear distorted with a long timescale. This 
behaviour may indicate a velocity variation owing to a slow 
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evolution of the stellar activity, rather than a long-period plan- 
etary companion. Both interpretations are possible however, be- 
cause the degree of correlation as shown in Figure 14 is not high 
(72%). 

When we fit the data in search for a Keplerian orbit, we find a 
circular orbit of a planet with a minimum mass of 0.6 Mj up and 
2.5 AU semi-major axis. It corresponds to a signal at a period 
of 1500 to 1900 d, or 4 to 5 years. This is compatible with the 
typical duration of an activity cycle, although on the short end 
(see e.g. Baliunas et al. (1995); Fares et al. (2009); Metcalfe et al. 
(2010)). 

3.5. HIP21934 

Twenty-three spectra of HIP 2 1934 were secured with HARPS 
from January 2004 to April 2010, with an average error of 
1.4 m/s. Low-amplitude velocity variations are detected, with 
clearly a long period variation and an additional short-period 
component (Figure 12). However, as for BD -114672, the full- 
width-at-half maximum of the cross-correlation function is cor- 
related to the activity index, with a coefficient of 89% (Figure 
14). This again indicates that the velocity variation may be 
caused by pure stellar activity fluctuations, with a slowly evolv- 
ing distortion of stellar lines over several years. In this case, the 
correlation is stronger, and the interpretation in terms of mag- 
netic cycle is favoured, although the star is moderately active 
during most of the cycle. The velocity variation of this cycle 
mimics the signature of a long -period (1100 days) planetary 
companion of minimum mass less than one Jupiter mass. This 
result is worse than expected in dedicated studies by Santos et al. 
(2010a), because the amplitude of the effect is higher - even 
though HIP 2 1934 is a K dwarf, as was considered in this study. 
It illustrates the importance of monitoring and considering ac- 
tivity tracers simultaneously to radial velocities for a valid inter- 
pretation of the results of planet search surveys, especially when 
low-velocity amplitudes and long orbital periods are considered. 
The length of the magnetic cycle, if confirmed as such, is 3-3.5 
years, again a rather short period compared with the sample of 
Baliunas et al. (1995). This would correspond to a ~ 7-yr mag- 
netic cycle, intermediate between the 2-yr cycle of tBoo (Fares 
et al. 2009) and 22-yr cycle of the Sun. 

4. Conclusion 

We presented the observational material secured with the high- 
precision radial-velocity spectrograph HARPS on ten stars 
in search for the signature of planetary systems. One star, 
HD 72659 , had a published planet, for which we could refine the 
mass determination and orbital elements. Figure 15 shows how 
these new planets hold on the total distribution of mass, period, 
and eccentricity of known exoplanets. 

After the first seven years of our survey, we aim at giving 
first synthetic views of the results by publishing our current 
detections of long period companions in the planetary mass 
regime. The radial-velocities time series reported in this work 
sometimes appear incomplete, and several different solutions 
cannot be distinguished, as illustrated for some cases by 
our x 2 distributions against related parameters P, K and e. 
Measurements are continuing however, and in a few years we 
will be able to better constrain the orbits and refine mass deter- 
mination of the depicted systems. The discovered companions 
are giant planets and one or two brown dwarfs, with orbital 
periods more than three years. In two cases, indications of 
velocity variations caused by stellar long-term activity evolution 
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Fig. 7. Radial-velocity measurements of HD 33473A obtained with 
HARPS against time (left) and phase (right). 
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Fig. 8. Radial- velocity measurements of HD 89839 obtained with 
HARPS against time. 



are shown, rather than planetary signatures. 

Finding long-period, massive planets in volume-limited 
surveys is of prime importance to get the complete picture 
of extrasolar system's architecture, despite the natural biases 
towards short-period planets of the radial-velocity method. 
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Fig. 9. Statistical tests on the RV time series of HD 89839 : (left) cor- Fig. 11. Statistical tests on the RV time series of HD 167677 : (left) cor- 
relation plot between period and eccentricity; (right) correlation plot relation plot between period and eccentricity; (right) correlation plot 



between velocity semi-amplitude and period. 



between velocity semi-amplitude and period. 
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Fig. 10. Radial-velocity measurements of HD 167677 obtained with 
HARPS against time. 



However, when searching for best-fit solutions of planetary 
orbits with periods larger than the time span of the observations, 
one meets the degeneracy of several types of solutions, with 
or without linear drifts, additional components in the system, 
or accounting for degeneracies in the fitted parameters. One 
expects long period planets to be part of systems, because in- 
ward migration should not have blown out the inner planets, and 
series of giant planets may exist as in the solar system. Orbital 
fitting for long-period planets is thus more uncertain than for 
shorter-period planets. The significance of periods, masses, and 



eccentricities of known radial-velocity long-distance planets 
should then be taken with some care, as the best-fit solution 
is likely to evolve when more data are available. One recent 
example of revisiting radial-velocity solutions of long-period 
planets of 47 UMa is illustrated by Gregory & Fischer (2010). 
The solutions given in the present work will also be more 
constrained with additional data, either to complete a period, 
or to confirm the presence of an additional component in the 
system. The data are still being collected with HARPS under 
the observing programme 085.C-0019 (P.I. Lo Curto) and will 
be presented and analysed in forthcoming papers. 

Direct imaging of the new planetary systems whose discov- 
ery is reported here will be possible, especially for the most 
nearby target HD 1 13538 . The angular separation presently ex- 
pected for their planets is 0.152 arcsec from its stars. 

As for other planets around giant stars, HD 1690 b is a 
massive gaseous planet, similar to the recently discovered 
HD 32518 b (Dollinger et al. 2009). There are now about 20 
known planetary systems to giant stars, all with semi-major 
axes exceeding 0.59 AU and minimum masses larger than 1 .5 
Mj llp . This narrower spread of parameters of planets around 
giant stars than around main-sequence stars may be explained 
by the high residual velocity caused by radial oscillations. 
With periods of hours to days and amplitude of several tens of 
m/s (Hatzes et al. 2007), these intrinsic variations of the star 
increase the difficulty of detecting the low-mass end of the 
planets. Regarding arguments such as the filling of the Roche 
lobe by the primary, the minimum period expected for a giant 
of l.O9M mass as HD 1690 is about 200 days (Mermilliod & 
Mayor 1996), which is below the observed 532 day orbit. 
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preliminary, these facts suggest that long period giant planets 
are not uncommon around low-metallicity stars, as shown by 
a dozen other examples (see http://exoplanet.eu/catalog.php). 
This possible trend will be confirmed by this continuing survey 
and others. 

Finally, let us stress that discoveries of giant planets in very 
long orbits is only beginning, because of the long baselines that 
are needed. When the star activity is not an issue, Jupiter twins 
will be soon collected by radial-velocity surveys, and such sys- 
tems will be good candidates when searching for Earth twins. 
The role of Jupiter-like planets as shields for the formation of 
inner terrestrial planets has been underlined decades ago (e.g. 
Wetherill (1991)). Finding the most massive planet in a system, 
on a distant orbit, may thus be a first step to discovering plane- 
tary systems similar to the solar system, and searching for Earth 
analogs. 



> 



■l.i.l.i.l.i.l. 



f : 

i . i . i . i . i . i . i . i . i 



3000 3200 3400 3600 3800 4000 4200 4400 4600 4800 5000 5200 5400 
JD - 2450000.0 [days] 

Fig. 12. Radial- velocity measurements against time of stars BD - 
1 14672 and HIP 21934 for which the long-term velocity variation is in- 
terpreted as caused by the magnetic activity cycle of the star. 



It may be a new fact that magnetic cycles of stars may have a 
signature on the radial-velocity behaviour over long timescales, 
with amplitudes similar to the gravitational perturbation of a gi- 
ant planet in a large-distance orbit. This was recently discarded 
by Santos et al. (2010a) for a limited sample of G and early K 
stars, and is further discussed by Lovis et al. (in prep). Two ex- 
amples are given here, of stars from the volume-limited sample 
followed-up by HARPS (HIP 21934 and BD -114672) showing 
the width of stellar lines affected by the long-term magnetic cy- 
cle rather than affected by a distant planetary companion. In this 
paper, we tentatively interpret the radial-velocity variations as 
caused by this long-term cycle, but investigations should be con- 
tinued in this field, from analyses of a larger sample and from 
simulations. 

One planetary system is presented around HD 113538 with 
a Saturne- Jupiter-like pair of planets at 0.7 and 2.4 AU respec- 
tively. Like the vast majority of planet pairs in systems, the most 
massive planet is the outer planet. This is possibly a detection 
bias towards the long periods, because the signal amplitude 
greatly decreases with star-planet distance. 

Interestingly, most of the planet host stars discussed here 
are metal-poor, contradicting the usual conclusion that stars 
with giant planets are usually more metal-rich than the average 
field dwarfs (e.g. Santos et al. (2004), Fischer & Valenti (2005). 
In a recent paper, Santos et al. (2010b) have also discovered 
three giant planets orbiting metal-poor stars from the HARPS 
survey, all in long period orbits. Similarly, it has been noted 
that the metallicity distribution of giant stars hosting planets 
does not favour metal-rich stars (Pasquini et al. 2007). Though 
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Table 1. Observed and inferred stellar parameters for the planet-hosting stars presented here. For stars with lowest masses, the relationships which 
derive the activity index and projected velocity are not calibrated, and the flux in the blue calcium lines is too faint for an estimate of log/?[j K . 



Parameter 


HD 1690 


BD -114672 


HD 25171 


HD 33473 A 


HD 72659 


HD 89839 


Sp 




K1III 


K7V 


F8V 


G3V 


G2V 


F7V 


V 


[mag] 


9.17 


10.02 


7.79 


6.71 


7.46 


7.64 


B-V 


[mag] 


1.354 


1.263 


0.554 


0.662 


0.612 


0.523 


n 


[mas] 


3.22 (1.43) 


36.65 (1.73) 


18.19(0.45) 


18.69 (0.49) 


20.07 (0.75) 


17.65 (0.55) 


d 


[pc] 


310(100) 


27.3 (1.3) 


55 (1.5) 


53.5 (1.5) 


49.8 (2) 


57 (2) 


My 


[mag] 


1.7 (0.7) 


7.84 (0.1) 


4.09 (0.07) 


3.07 (0.06) 


3.97 (0.1) 


3.86 (0.07) 


B.C. 


[mag] 


-0.75 


-0.61 


-0.03 


-0.09 


-0.06 


-0.015 


L 


[£o] 


33.1 


0.1 


1.89 


5.11 


2.17 


2.30 


Teg 


[K] 


4393 (85) 


4594 (134) 


6160 (65) 


5740 (62) 


5926 (62) 


6314 (65) 


logg 


[cgs] 


2.12(0.17) 


4.34 (0.31) 


4.43 (0.1) 


3.97(0.1) 


4.24(0.1) 


4.49 (0.12) 


[Fe/H] 


[dex] 


-0.32 (0.06) 


-0.47 (0.08) 


-0.11(0.044) 


-0.13 (0.04) 


-0.02 (0.098) 


+0.04 (0.043) 


M, 


[M ] 


1.09 (0.15) 


0.7 (0.1) 


1.09 (0.03) 


1.23 (0.03) 


1.06 (0.04) 


1.21 (0.03) 


v sin i 


[km s- 1 ] 


3.5 




1.0 


2.6 


1.9 


4.1 


S MW 




0.11 (0.04) 


0.75 (0.06) 


0.15 (0.01) 


0.15 (0.01) 


0.16(0.01) 


0.15 (0.015) 


l°g*HK 








-4.99 


-5.12 


-4.97 


-4.97 


Age 


[Gy] 


6.7 (3.2) 


4. (3.7) 


4.0(1.6) 


4.4 (0.3) 


6.5 (1.5) 


1.2 (0.9) 


R, 


[R ] 


16.7 (3.6) 


0.7 


1.18(0.04) 


2.23 (0.09) 


1.36 (0.06) 


1.25 (0.04) 



Table 2. Observed and inferred stellar parameters for the planet-hosting stars presented in this paper, (cont'd) 



Parameter 


HD 113538 


HD 167677 


HD 217786 


HIP 21934 


Sp 




K9V 


G5V 


F8V 


K3 


V 


[mag] 


9.02 


7.90 


7.80 


9.9 


B-V 


[mag] 


1.362 


0.705 


0.578 


1.197 


n 


[mas] 


63.03 (1.36) 


17.29(1.04) 


18.23 (0.72) 


28.77 (0.97) 


d 


[pc] 


15.8 (0.4) 


58(3) 


54.8 (2) 


34.8 (1.2) 


M v 


[mag] 


8.03 (0.06) 


4.08 (0.1) 


4.11 (0.07) 


7.19(0.08) 


B.C. 


[mag] 


-0.66 


-0.135 


-0.05 


-0.62 


L 


lL a ] 


0.09 


2.09 


1.89 


0.19 


T s ff 


[K] 


4685 (155) 


5474 (65) 


5966 (65) 


4674 (172) 


logg 


[cgs] 


4.28 (0.37) 


4.43 (0.1) 


4.35 (0.11) 


4.2 (0.41) 


[Fe/H] 


[dex] 


-0.17 (0.074) 


-0.29 (0.043) 


-0.135 (0.043) 


+0.025 (0.12) 


M, 


[M G ] 


0.7 (0.1) 


0.96 (0.02) 


1.02(0.03) 


0.72 (0.02) 


v sin i 


[km s-»] 




< 1 


1.4 




S MW 




1.05 (0.08) 


0.17(0.015) 


0.15 (0.01) 


0.57 (0.06) 


l0g^ K 






-4.99 


-4.86 


-4.86 


Age 


[Gy] 


1.3(1.3) 


11(1) 


7.6 


3.9 (3.8) 


R, 


[Ro] 


0.65 (0.1) 


1.7 (0.1) 


1.27 (0.06) 


0.7 (0.1) 
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Table 3. Orbital and physical parameters for the planets presented in this paper. T is the epoch of periastron. tr(O-C) is the residual noise after 
orbital fitting of the combined set of measurements, linear is the slope of the observed drift with time. 



Parameter 


HD 1690 b 


HD 25171 b 


HD 72659 b 


HD 217786 b 


P 


[days] 


533 (1.7) 


1845 (167) 


3658 (32) 


1319(4) 


T 


[JD-2400000] 


54449 (5) 


55301 (449) 


55351 (59) 


56053 (7) 


e 




0.64 (0.04) 


0.08 (0.06) 


0.22 (0.03) 


0.40 (0.05) 


linear 


[m/s/yr] 


-7.2 (0.4) 


- 


- 


- 


J 


[km s" 1 ] 


18.23 (0.03) 


43.6 (0.01) 


-18.15 (0.001) 


9.98 (0.006) 


OJ 


[deg] 


122 (8) 


96 (89) 


-99.3 (6.5) 


100 (2) 


K 


[m s" 1 ] 


190 (29) 


15.0 (3.6) 


41.0(1.3) 


261 (20) 


a\ sin i 


[10~ 3 AU] 


7.1(1) 


2.5 (0.3) 


13.4 (0.4) 


28.8(1.5) 


f(m) 


[10- 9 M o ] 


181 (84) 


0.66 (1.2) 


24.2 (2.2) 


1857 (307) 




[M Jup ] 


6.1 (0.9) 


0.95 (0.10) 


3.15 (0.14) 


13.0 (0.8) 


a 


[AU] 


1.3 (0.02) 


3.02 (0.16) 


4.74 (0.08) 


2.38 (0.04) 


N meas 




42 


24 


50a 


15 


S pan 


[days] 


2547 


2314 


4515 


1538 


cr (O-C) 


[m s- 1 ] 


29 


3.4 


4.2 


2.5 


X re d 




1.6 


1.9 


3.7 


2.1 



Table 4. Orbital and physical parameters for HD 1 13538 b and c 1 . 



Parameter 


HD 113538 b 


HD 113538 c 


P 


[days] 


263.3 (2.3) 


1657 (48) 


T 


[JD-2400000] 


54950 (8.5) 


56110(515) 


e 




0.61 (0.11) 


0.32 (0.06) 


OJ 


[deg] 


74 (20) 


-90 (21) 


K 


[m s->] 


15.5 (12) 


16.5 (1.4) 


a\ sin i 


[10~ 3 AU] 


0.26 (0.07) 


2.36(0.15) 


f{m) 


[1O" 9 M ] 


0.048 (0.10) 


0.65 (0.14) 


m 2 sin i 


[M Jup ] 


0.27 (0.08) 


0.71 (0.06) 


a 


[AU] 


0.71 (0.01) 


2.43 (0.06) 


J 


[kms- 1 ] 


39.23 (0.012) 




Ameas 




29 




S pan 


[days] 


2268 




cr (O-C) 


[m s- 1 ] 


1.8 




V 2 
A red 




1.5 





Table 5. Orbital and physical parameters for the companions with incomplete orbits. 



Parameter 


HD 33473 A b 


HD 89839 b 


HD 167677 b 


P 


[days] 


2798 (60) 


6601(-3570,+4141) 


1814(100) 


T 


[JD-2400000] 


55053 (93) 


60067 (-4954,+4425)) 


55656 (123) 


e 




0.07 (0.01) 


0.32 (0.2) 


0.17(0.07) 


linear 


[m/s/yr] 


-78.9 (0.07) 






J 


[kms" 1 ] 


44.2 (0.001) 


31.7 (0.002) 


-57.24 (0.01) 


OJ 


[deg] 


141 (15) 


-22 (144) 


-80 (24) 


K 


[m s~'] 


91 (1.1) 


45.5 (2.3) 


23.7 (3) 


a\ sin; 


[lO- 3 AU] 


23.4 (0.7) 


22 (17) 


3.9 (0.4) 


f(m) 


[lO-'Mo] 


218 (12) 


38 (10) 


2.4 (2.3) 


m 2 sin ( 


[M Jup ] 


7.2 (0.3) 


3.9 (0.4) 


1.36 (0.12) 


a 


[AU] 


4.17 (0.09) 


6.8 (-2.4,-1-3.3) 


2.9 (0.12) 


^meas 




26 


39 


26 


5 pan 


[days] 


2465 


2297 


1977 


cr (O-C) 


[m s- 1 ] 


2.4 


3.0 


2.6 


v 2 

A red 




3.5 


1.9 


3.1 
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Fig. 13. Bisector slope variations with respect to radial velocity for all data sets. 
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Fig. 14. Variations of the activity tracer, logR' HK or Mount-Wilson 5 index, are correlated to the full-width-at-half maximum of the cross-correlation 
function of stars BD -1 14672 (top) and HIP 21934 (bottom). This indicates that the long-term velocity variation is probably caused by the magnetic 
cycle of the star, rather than a long-period planetary companion. The correlation is more significant for star HIP 21934 . 
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Fig. 15. Mass-period and period-eccentricity diagrams. Open circles show the full sample of radial-velocity planets, as of July 2010, and the filled 
circled depict the planets described here. Green filled circles related by a line show the location of planet HD 33473A b considering both solutions 
that are equivalent (see text); only one of those dots is a planet. 



